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Abstract. The lattice QCD studies indicate that the critical temperature Tc ~ 260 — 280 MeV of the 
deconfinement phase transition in quenched QCD is considerably smaller than the lowest-lying glueball 
mass mc — 1500 — 1700 MeV, i.e., Tc <^ mo. As a consequence of this large difference, the thermal 
excitation of the glueball in the confinement phase is strongly suppressed by the statistical factor as 
g-™-G/Tc ^ 0.00207 even near T ~ Tc. We consider its physical implication, and argue the abnormal feature 
of the deconfinement phase transition in quenched QCD from the statistical viewpoint. To appreciate this, 
we demonstrate a statistical argument of the QCD phase transition using the recent lattice QCD data. 
From the phenomenological relation among Tc and the glueball mass, the deconfinement transition is found 
to take place in quenched QCD before a reasonable amount of glueballs is thermally excited. In this way, 
quenched QCD reveals a question "what is the trigger of the deconfinement phase transition ?" 



PACS. 12.38.Mh Quark-gluon plasma - 12.38.Gc Lattice QCD calculations 
12.39.Mk Glueball and nonstandard multi-quark/gluon states 



12.39.Ba Bag model 



1 Introduction 

The quark-gluon-plasma (QGP) is one of the most in- 
teresting targets in the finite-temperature quark-hadron 
physics I] . Currently, the QGP creation experiment is be- 
ing performed in the relativistic heavy ion collider (RHIC) 
project at the Brookhaven National Laboratory (BNL), 
and much progress in understanding the finite-temperature 
QCD is desired. Historically, the instability of the hadron 
phase was first argued by Hagedorn |2| before the dis- 
covery of QCD. He pointed out the possibility of a phase 
transition at finite temperature, based on the string or the 
flux-tube picture of hadrons|2I21- After QCD was estab- 
lished as the fundamental theory of the strong interaction, 
this phase transition was recognized as the deconfinement 
phase transition to the QGP phase, where quarks and 
gluons are liberated with the restored chiral symmetry. 
The QCD phase transition has been studied using var- 
ious QCD-motivated effective models such as the linear 
a model@], the Nambu-Jona-Lasinio model 5 , the dual 
Ginzburg-Landau theory|H] and so on. 

In order to study nonperturbative features of the QCD 
phase transition, the lattice QCD Monte Carlo calcula- 
tion serves as a powerful tool directly based on QCD. It 
has been already extensively used to study the nature of 
the QCD phase transition. At the quenched level, SU(3) 
lattice QCD indicates the existence of the deconfinement 



phase transition of a weak first order at Tc ~ 260 — 280 
MeV 7 . On the other hand, in the presence of dynami- 
cal quarks, it indicates the chiral phase transition at Tc = 
173(3) MeV for Nj = 2 and Tc = 154(8) MeV for iV/ = 3 
in the chiral limit 8 . 

For the comparison with the experimental data in the 
real world, it would be desirable to investigate full QCD 
with dynamical quarks. However, there are a number of 
underlying nonperturbative features, which are shared in 
common by both full QCD and quenched QCD such as 
color confinement and instanton phenomena. In order to 
understand such nonperturbative features of QCD, quenched 
QCD provides us with an idealized environment to fo- 
cus on the essence of the problem itself without involving 
inessential technical complexities. 

The aim of this paper is to point out the abnormal fea- 
ture of the large difference between the lowest-lying glue- 
ball mass mo — 1500 — 1700 MeV and the critical tem- 
perature Tc = 260 - 280 MeV in quenched QCD. In fact, 
due to this large difference, the thermal exciation of the 
glueball is strongly suppressed by a small statistical factor 
g-mc/T, ^ 0.00207 even near T ~ Tc. In order to appre- 
ciate this point, we demonstrate a statistical argument on 
the QCD phase transition using the closed packing model 
with the bag-model picture of hadrons. The closed pack- 
ing model has been often used for the phenomenological 
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understanding of the full QCD phase transition near Tc 
from below. In the confinement phase, only color-singlet 
states such as hadrons can contribute to the Boltzmann 
sum in the partition function, and therefore it is desirable 
to understand the QCD phase transition in terms of the 
hadronic degrees of freedom below Tc. 

The contents are organized as follows. In Sect. [3 we 
give a brief review of the closed packing model, i.e., the 
statistical treatment of the QCD phase transition with the 
bag-model picture of hadrons. We derive a phenomenolog- 
ical relation among the hadron masses, the hadron sizes 
and the critical temperature. In Sect. 13 we apply the sta- 
tistical approach to full QCD and compare it with the 
recent lattice data of full QCD. In Sect. 0] we apply the 
statistical approach to quenched QCD, and demonstrate 
an abnormal feature of the quenched QCD phase tran- 
sition. In Sect. 13 after the summary of the results, we 
attempt to clarify the essence of this problem, and discuss 
an abnormal nature of the deconfinement phase transition 
in quenched QCD. 
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Fig. 1. The characteristic function /(m) on the thermal boson 
in Eq. Q plotted against m = m/T. 



2 A review of the statistical approach 

In this section, we give a brief review of the closed pack- 
ing model, i.e., the statistical treatment of the QCD phase 
transition based on the bag-model picture of hadrons. It 
provides us with a simple but useful insight into the QCD 
phase transition from below Tc in terms of hadronic prop- 
erties. We derive a phenomenological relation between the 
critical temperature Tc and the properties of the hadrons, 
i.e., the mass and the size. 

In the bag-model picture, quarks and gluons are as- 
sumed to be confined inside a spherical bag. Color con- 
finement is simply taken into account through the bag- 
like intrinsic structure of hadrons jH]. At low temperature, 
only a small number of such bags are thermally excited, 
and the thermodynamic properties of the system are de- 
scribed in terms of these spatially isolated bags. With the 
increasing temperature, the number of the thermally ex- 
cited bags increases. Gradually, these bags begin to over- 
lap one another, and they finally cover the whole space 
region at a critical temperature Tc in this picture. Above 
Tc, as a result of overlapping bags, the whole space is filled 
with liberated quarks and gluons, and the thermodynamic 
properties of the system is now governed by quarks and 
gluons. In this way, the QCD phase transition is described 
in terms of the overlaps of the thermally excited bags. 

For the quantitative argument, we define the spatial 
occupation ratio ry (T) at temperature T to be the ratio of 
the total volume of the spatial regions inside the thermally 
excited bags to the volume V of the whole space region. In 
the closed packing picture of the QCD phase transition, 
rv{T) plays the key role, which is estimated as 



rv{T) 



V 



(1) 



(Pk 



where Nn{T), A„, m„ and i?„ are the number at temper- 
ature T, the degeneracy, the mass and the bag radius of 
the rt-th elementary excitation, respectively. Here, /(m) 
is defined by 



47r f d^k 



1 



1 
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and its functional form is plotted against m = m/T in 
Fig- n Note that f{m/T) is a characteristic function to 
describe the thermal contribution of the boson with the 
mass TO at the temperature T JUj. For m ^ T, f{m/T) 
decreases exponentially with m/T, and the thermal con- 
tribution is expected to become negligible. 

In the closed packing picture, the phase transition is 
assumed to takes place, when the thermally excited bags 
almost cover the whole space region. Hence, the critical 
temperature Tc is estimated by solving 



rv{Tc) = 1. 



(3) 



(27r)3 y/jnl+k-^/T 



In the closed packing picture, an essential role is played 
by the color confinement, which is a peculiar phenomenon 
in QCD making the underlying quark and gluon struc- 
ture hidden inside the hadron. This feature is quite dif- 
ferent from atomic system. (For instance, the closed pack- 
ing picture of phase transition cannot be applied to the 
Coulombic system such as the ionization transition of an 
atomic gas, where the bag-model description of the bound- 
state is inappropriate.) The readers might feel that the 
closed packing argument is too naive for the complicated 
QCD phase transition. In particular, the ideal-gas statisti- 
cal treatment should be applied only to the dilute system 
in a strict sense. Nevertheless, the closed packing approach 
provides a simple and physical insight into the complicated 
QCD phase transition in terms of the thermal excitations 
of hadrons. Hence, it would be natural to attempt under- 
standing the QCD phase transition in the closed packing 
picture as a first step. 
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Fig. 2. The critical temperature Tc plotted against the pion 
mass niTr in the Nf = 2 case. The dashed curve denotes the 
result retaining only the contribution from pions. The solid 
curve denotes the result including also contributions from p 
and uj mesons. The triangle denotes the lattice data taken 
from Ref.(Ilj obtained with standard (std) fermion action. The 
square denotes the lattice data taken from Ref.[S| obtained 
with the improved staggered (p4) fermion action. 
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Fig. 3. The critical temperature Tc plotted against the pseudo- 
scalar meson mass mps in the SU(3)/ symmetric case. The 
dashed curve denotes the result retaining only the contribu- 
tion of the pseudo-scalar mesons. The dot-dashed and the solid 
curves denote the results including also the contribution of the 
flavor-octet vector mesons with i? = 1 fm and R = 0.75 fm, 
respectively. The circle denotes the lattice data taken from 
Ref.|H] obtained with the improved staggered (p4) fermion ac- 
tion. 



3 A statistical approach to the full-QCD 
phase transition 

First, we apply the statistical approach with the closed 
packing model to the full QCD phase transition, and com- 
pare the results with the recent lattice QCD data. In full 
QCD, the lightest physical excitation is the pion, and all 
the other hadrons are rather heavy as m ^ m^, Tc. In fact, 
the pion is considered to play the key role in describing 
the thermodynamic properties of full QCD below Tc from 
the viewpoint of the statistical physics. Hence, in most 
cases, only the pionic degrees of freedom are taken into 
account in the hadron phase in the argument of the full- 
QCD phase transition. By using the isospin degeneracy 

= 3, the mass = 140 MeV and the radius Rt^ — 1 
fm, we solve Eq. © with rv{T) = 3RlT^f{m^/T) to 
estimate the critical temperature as Tc ~ 183 MeV. Con- 
sidering its closeness to the full lattice QCD result with 
Nf = 2, i.e., Tc ~ 170 MeV, the statistical approach to 
the full-QCD phase transition seems to be rather good. 

We now consider the m7r-dependence of the critical 
temperature Tc. Note that, in the actual lattice QCD cal- 
culations, the pion mass is taken to be still rather heavy 
as TTT-TT > 400 MeV for the technical reasons. From these 
data, the critical temperature Tc in the chiral limit is ob- 
tained using the chiral extrapolation. In Ref.|H], the au- 
thors parametrized the m,r-dependence of Tc in the full 
lattice QCD with TV/ = 2 as 



^) M - 0.40(1) + 0.039(4) (^^^ , 



(4) 



where a denotes the string tension. Strictly speaking, the 
phase transition becomes just a cross-over for intermedi- 
ate values of m-^. Hence, in Ref.jHI, the pseudo critical 



temperature is adopted as Tc, which is determined from 
the peak positions of the susceptibilities of the Polyakov 
loop and so on. 

For the Nf = 2 case. Fig. |21 shows the estimate of Tc 
plotted against ttItt based on Eq. |5J) in the closed packing 
model. The dashed curve denotes Tc retaining only the 
contribution from pions with R^^ ~ 1 fm. The solid curve 
denotes the result in the case including also the contri- 
butions from the low- lying vector mesons, such as p and 
u!, which are the next lightest particles in Nf — 2 full 



QCD. Here, we have used Ap = 3 x 3 = 9, mp 



770 



MeV, Rp ~ 1 fm, A„ = 3, m„ = 783 MeV, i?^ ~ 1 fm as 
inputs, which are treated as mjr-independent constants. 
These vector mesons give an additional contribution to 
rv{T) as 6rv{T) = m^T^ f{mp/T) + SR^T^ fim^/T). 
The triangle and the square in Fig. |21 denote the lattice 
data taken from Refs. pilllj . The m7r-dependence of the 
critical temperature Tc is improved after the inclusion of 
the low-lying vector mesons, i.e., p and uj. In spite of the 
naively mindedness of the closed packing model, this sta- 
tistical analysis seems to reproduce the lattice data, which 
may be surprising. 

We next consider the idealized SU(3) / symmetric case. 
In this case, the pseudo-scalar (PS) octet mesons, such as 
pions, kaons and i]s, are the lightest and possess the same 
mass mps in common with the degeneracy as Aps = 8. We 
plot, in Fig. |31 the critical temperature Tc against mps- 
The dashed curve denotes Tc retaining only the contri- 
bution from PS-mesons with i?ps — 1 fm, which leads 
to rv{T) = 8RlsT^f{mps/T). The dot-dashed curve in 
Fig. 121 denotes the results including also the contributions 
from the octet vector mesons with Ay = 8x3 = 24, 
my ~ 770 MeV, i?v — 1 fni. (Inclusion of the flavor- 
singlet vector meson does not change the result so much.) 
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These vector mesons give an additional contribution to 
rv{T) as Srv{T) = 2m%T^ f{mY /T). The circle in Fig.El 
denotes the lattice data taken from Ref.jS]. We see that 
both the dashed and the dot-dashed curves are roughly 
consistent with the lattice QCD results. We note that the 
small deviation almost disappears by slightly adjusting 
the bag size as i?ps — Rv ~ 0.75 fm, as is shown in Fig. 13 
with the solid curve. 

Thus, the simple statistical approach with the closed 
packing picture is seen to reproduce the recent lattice 
QCD data on the critical temperature Tc and its mps- 
dependence in the full QCD phase transition both for 
Nf = 2 and Nf = 3, which would suggest that some of 
essential natures of the full QCD phase transition could 
be governed by low- lying hadrons. Considering its naivete 
and its simply minded nature, this coincidence may be sur- 
prising. It would be interesting to refine this approach by 
including the interaction among hadrons. However, we em- 
phasize that the application of the closed packing model 
to full QCD is not the final aim of this paper. Rather, our 
aim is to demonstrate the mysterious mismatch between 
the critical temperature and the mass of the elementary 
excitation mode in the quenched QCD phase transition in 
the next section. 



4 A mystery in the quenched QCD phase 
transition 

In this section, we come to the main point of the pa- 
per, i.e., the crucial mismatch of the critical temperature 
Tc and the low-lying elementary-excitation mode in the 
quenched QCD phase transition. Focusing on the large 
difference between the lowest-lying glueball mass toq — 
1.5 - 1.7 GeV and Tc ~ 280 MeV in quenched QCD, we 
consider its physical implications. To make the argument 
more quantitative, we combine the recent lattice QCD 
data with the statistical approach, and attempt to make 
an estimate of the critical temperature Tc- 

In quenched QCD, due to the color confinement, only 
the color-singlet modes can appear as physical excitations. 
Since all of them consist of color singlet combinations of 
gluons, they are called as glueballs. The mass spectrum 
of the glueballs is known through the quenched lattice 
QCD calculations I12||13|IT1I . The lightest physical excita- 
tion is the 0+"'" glueball with toq = 1.5 — 1.7 GeV. Being 
the lightest physical excitation, the lowest 0+"'' glueball is 
expected to play the key role in describing the thermo- 
dynamic properties of quenched QCD in the confinement 
phase. Hence, we first take into account only the 0+"'" glue- 
ball. As the mass mG(s) and the size -Rg(S) of the lowest 
0++ glueball, we adopt the recent lattice QCD results as 
mG(s) = 1730 MeV [El and i?G(s) = 0.4 fm [TaTCj . We 
use these values together with the degeneracy Ag(s) = 1 
as inputs. Then the spatial occupation ratio is given by 
ry(T) — i?G(S) /("^G(S)/2^)i ^-iid the critical tempera- 
ture is estimated as Tc ~ 827 MeV from Eq. (jSJ. This es- 
timate is too much larger than the quenched lattice QCD 
result as Tc = 0.629(3) ~ 280 MeV in Ref. 7 with 
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Fig. 4. The critical temperature Tc plotted against the lowest 
0^''' glueball mass mQ(^s) in the statistical approach. The solid 
line denotes Tc for the glueball size -Rg(S) = 0.4 fm, which 
is taken from the recent lattice QCD [T?llfci| . The cross (x) 
indicates the quenched lattice QCD results, i.e., Tc = 280 MeV 
and mG(s) — 1730 MeV. We add the cases with -Rg(S) = 0.8 
and 1.2 fm by the dashed and dotted lines, respectively. 



y/a — 450 MeV. (In other words, only a tiny fraction of 
the space region can be covered by the thermally excited 
bags as ry(r) = 0.0021 at T = 280 MeV.) This large 
discrepancy would be a serious problem in the quenched 
QCD phase transition. 

To seek the solution, we examine the possibilities of 
the polemass reduction and the thermal swelling of the 
glueball at T ~ Tc, which are suggested in Ref.|jH|. We 
first consider the possibility of the polemass reduction at 
T ~ Tc. To this end, we show in Fig. 01 the critical tem- 
perature Tc plotted against the lightest 0++ glueball mass 
mG(s) in tli6 closed packing model for the glueball size 
i?G(S) = 0.4 fm with the solid curve. The cross (x) indi- 
cates the quenched lattice QCD results, Tc — 280 MeV 
and Tog — 1730 MeV. We see that, for the problem to 
be settled, the polemass reduction must be as significant 
as tog(s)(7c) ;$ 500 MeV. However, in the recent lattice 
QCD calculations [TMT?)] . it has been reported that the 
thermal 0++ glueball persists to hold a rather large pole- 
mass as tog(s)(2^ — Tc) ~ 1250 MeV. Hence, we examine 
another possibility, i.e., the thermal swelling. In Fig. ^ 
we also include Tc corresponding to different sizes of the 
glueball as i?G(S) = 0.8, 1.2 fm, which are denoted by the 
dashed and dotted curves, respectively. We see that to re- 
produce Tc ~ 280 MeV, the thermal swelling of the glue- 
ball must be quite significant at T ~ Tc. The explicit 
calculation leads to an abnormally large size i?G(S) — 3.1 
fm. However, such a drastic thermal swelling of the low- 
est scalar glueball was rejected by the recent lattice QCD 
studies [T5til5) , which states that the thermal glueball size 
is almost unchanged even near Tc. 

We seek for another possibility by including the con- 
tribution of the excited-state glueballs. In addition to the 
lowest 0++ glueball with tog(s) — 1730 MeV and i?G(S) — 
0.4 fm, we consider the thermal contribution from the 
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lowest 2++ glueball, which is the next lightest hadron in 
quenched QCD. We take Ag(t) = 5, tog(t) - 2400 MeV 
|12|. i?G(T) — 1 fin- (Here, we assume it to have a typical 
hadron size.) Then, the spatial occupation ratio receives 
a correction as Srv{T) = bR^^rj^s^T^ f {mQ(^T^ /T) , and the 

correction amounts to (5ry(T) = 0.0223 at T = 280 MeV. 
The resulting critical temperature is given as Tc ~ 432 
MeV, which is still too large. We note that the realistic 
glueball size would be more compact. However, if so, its 
contribution becomes more negligible. Besides these two 
low-lying glueballs, the following excited states are pre- 
dicted in Ref.[I2] as 0-+(2590), 0*++(2670), l+-(2940), 
2-+(3100), 3+-(3550), 0*-+(3640), 3++(3690), f — (3850), 
2*-+(3890), 2— (3930), 3— (4130), 2+-(4140), 0+^(4740). 
We include all the contribution from these excited states, 
assuming the unknown glueball size as a typical hadron 
size, i.e., Rq ~ 1 fm. The correction by these excited states 
amounts to only 5rv{T) = 0.0113 at T = 280 MeV, and 
the resulting critical temperature is estimated as — 395 
MeV, which is again too large. Note that we have used 
Rq ~ 1 fm for the radii of the excited glueballs, since 
we adopt it as a typical hadron size. However, even if we 
adopt Rq ~ 2 fm, which may be considered well beyond 
the upper bound for a single hadron radius, the total con- 
tribution to rviT) is 5rv{T) = 0.269 at T = 280 MeV, 
which is still insufficient. 

One may argue that our statistical argument can be 
improved by considering the self-interaction of the glue- 
ball and the Lorentz contraction effects. However, unlike 
in full QCD, the glueball system at T ~ Tc is so dilute 
that ideal-gas statistical treatment can be applied. (We 
will come back to this point in the summary.) On the 
other hand, due to the large difference of and the glue- 
ball masses, the effect of the Lorentz contraction is neg- 
ligible. Furthermore, the Lorentz contraction effect drives 
the discrepancy in the unwanted direction. 

To recapitulate this section, we have observed that the 
statistical approach leads to a terrible overestimate of Tc 
in quenched QCD even after the so many improvement 
have been attempted. The direct cause of this failure is 
the extremely small statistical factor e~™c:(s)/T'c ^ 0.0021, 
which strongly suppresses the excitations of the glueballs. 
As a consequence, only a small fraction of the space region 
can be covered by the thermally excited bags of glueballs 
even at T = 280 MeV, and Tc becomes abnormally large. 

5 Summary and discussions — What is the 
trigger or the driving force of the QCD phase 
transition ? 

We have considered the large difference between the glue- 
ball mass ruQ ~ 1500 — 1700 MeV and the critical tem- 
perature Tc ~ 260 - 280 MeV in the quenched QCD. As 
a consequence of this large difference, the thermal exci- 
tation of a single glueball is suppressed by a strong sta- 
tistical factor ^-"^gITc ^ 0.00207. We have considered its 
physical implications and argue the abnormal nature of 
the quenched QCD phase transition. To appreciate how 



abnormal it is, we have used the statistical argument with 
the bag-model picture of hadrons, i.e., the closed pack- 
ing model. We have derived a phenomenological relation 
among the critical temperature Tc, the mass and the size 
of the low- lying hadrons. We have demonstrated that with 
slight modifications of the model parameters, the closed 
packing model can reproduce the mTr-dependence of Tc ob- 
tained by full lattice QCD for TV/ = 2 and 3 cases, suggest- 
ing that some of the essential ingredients of the full QCD 
phase transition is governed by the lowlying hadrons. 

Unlike full QCD, we have found that the statistical 
approach terribly overestimates the critical temperature 
as Tc — 827 MeV in quenched QCD, and that only a tiny 
fraction of the space region is covered by the thermally 
excited bags of glueballs at T ~ 280MeV. We have consid- 
ered the possibility of the thermal swelling of the glueball 
size and the polemass reduction of the glueball near the 
critical temperature. However, both of these two have not 
provided us with the solution on the large discrepancy. 
Even though one includes all the contributions from the 
15 low- lying glueballs up to 5 GeV predicted in quenched 
lattice QCD 12 , the discrepancy remains to be still large. 
In other words, the number of the thermally-excited glue- 
balls is still too small at T ~ 280MeV, even after so many 
glueball excited states are taken into account. The direct 
origin of this discrepancy is the strong suppression of the 
thermal excitation of glueballs due to the extremely small 
statistical factor as e-™G(s)/T ^ 0.00207 at T = 280 MeV 
even for the lightest glueball, which leads to the insuffi- 
cient amount of the covered space by the thermally excited 
bags of glueballs. We remark that the crucial role is played 
by the smallness of the statistical factor e~™"2(s)/-^, which 
has rather general nature. Hence, through this failure, the 
closed packing argument suggests a general tendency that, 
if based on the lowlying excitation modes, the natural 
critical temperature Tc tends to become much larger in 
quenched QCD. Actually, such a tendency is also found 
in the field-theoretical model such as the dual Ginzburg- 
Landau theory [HI, which is an effective theory of color 
confinement based on the dual Meissner effect. It is re- 
markable that the essential origin of this large deviation 
of Tc can be understood in our naive statistical approach 
in a simplified and idealized manner. 

Although several arguments given so far have been 
based on the bag-model picture, this problem itself has 
a quite general nature, which can go beyond the reliabil- 
ity of the model framework. Finally, we reformulate this 
problem in a model-independent general manner by con- 
sidering the inter-particle distance I instead of the bag 
size R. Note that the reliability of the statistical argument 
becomes improved in the dilute glueball gas limit. Now, 
taking into account all the low-lying 15 glueball modes up 
to 5 GeV predicted in quenched lattice QCD we cal- 
culate the inter-particle distance I of the glueballs based 
on only the statistical argument. At T = 280 MeV, the 
inter-particle distance is estimated as / ~ 5 fm. It fol- 
lows that the deconfinement phase transition takes place 
at T ~ 280MeV, where the glueball density is rather small 
as p = l/{^(2.5fm)3} ~ l/(4.0fm)3. Since the theo- 
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retical estimates of the glueball size are rather small as 
Rg{s) < 0.4 fm jl5P6"17', the glueball system can be re- 
garded to be dilute enough at this density just below Tc. 
Furthermore, since the long-range interaction among glue- 
balls is mediated by the virtual one-glueball exchange pro- 
cess in quenched QCD, the interactions among glueballs 
are exponentially suppressed beyond its Compton length 
l/mG(s) = 0.112 fm. Therefore, one cannot expect the 
strong long-range interaction acting among the spatially- 
separated thermal glueballs, and the dilute gas treatment 
of the glueballs is considered to be valid at least in the 
confinement phase in quenched QCD. 

Now, it is quite difficult to imagine how such a too rare 
excitation of thermal glueballs can lead to the phase tran- 
sition. We arc thus arrive at the mystery. What is really 
the trigger or the driving force of the deconfinement phase 
transition in quenched QCD ? Actually, also in full QCD, 
the trigger of the QCD phase transition would be the in- 
teresting quenstion, but not so many have been known 
yet. In this sense, the careful consideration on this prob- 
lem may provide us with a key to find out the trigger of 
the QCD phase transition. 

H. S. is supported by Grant for Scientific Research (No. 12640274) 
from Ministry of Education, Culture, Science and Technology, 
Japan. 
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